A semiclassical theory of the Anderson transition 
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We study analytically the metal-insulator transition in a disordered conductor by combining the 
self-consistent theory of localization with the one parameter scaling theory. We provide explicit 
expressions of the critical exponents and the critical disorder as a function of the spatial dimen- 
sionality d. The critical exponent v controlling the divergence of the localization length at the 
transition is found to be v = | + -jr^ thus confirming that the upper critical dimension is infinity. 
Level statistics are investigated in detail. We show that the two level correlation function decays 
exponentially and the number variance is linear with a slope which is an increasing function of the 
spatial dimensionality. Our analytical findings are in agreement with previous numerical results. 

PACS numbers: 72.15.Rn, 71.30.+h, 05.45. Df, 05.40.-a 
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The recent experimental realization of disorder in ultra 
cold atoms 0] together with the rapid progress in numer- 
ical calculations [1, d, 0, H, 0, 0] has revived the interest 
in the metal-insulator transition (MIT) 0]. The statis- 
tical analysis of the spectrum and eigenvectors plays a 
central role in the identification and characterization of 
a MIT. Typical signatures of a MIT include: a) multi- 
fractal 0,0] eigenstates (for a review see 0,E(1), namely, 
the scaling of V q = j d d r\ijj(r)\ 2q oc L-^fa- 1 ) with re- 
spect to the sample size L is anomalous with D q < d 
a set of exponents describing the transition, b) scale in- 
variance [3( of the spectral correlations, c) level repul- 
sion of neighboring eigenvalues as in a disordered metal, 
d) linear number variance like in a disordered insulator, 
= ((Ni - (N £ )) 2 ) ~ X t (N t is the number of eigen- 
values in an interval of length < > 1 and (...) stands 
for ensemble average) but with a slope x < !• I n 3d, 
X ~ 0.27 for higher dimensions see Ref.0, 0j- Level 
statistics with these features are usually referred to as 
critical statistics [3. 11. 



Theoretical progress has been much slower in recent 
times. With the exception of the case of a Cayley tree 
[IH geometry we are indeed far from a quantitative an- 
alytical theory of the MIT. Below we briefly review the 
main analytical approaches to the MIT problem. In the 
original Anderson's paper the critical disorder at which 
the MIT occurs was estimated by looking at the limits 
of applicability of a locator expansion 0, In this 

approach 0] the metal-insulator transition is induced by 
increasing the hopping amplitude of an initially localized 
particle. A more refined self-consistent fl5| condition, 
still within the locator formalism, provided with a sim- 
ilar answer. This method is only exact in the case of a 
Cayley tree but it is believed to be accurate in the d — > oo 
limit. In both cases a MIT in 3d is predicted correctly. 
However the estimated critical disorder is considerably 
smaller than the one found in numerical simulations 0]. 
This disagreement persists in higher dimensions 0]. 

The one parameter scaling theory [l(| (OPT) pro- 
vides a completely different approach to the MIT. A key 
concept in this theory is the dimensionless conductance 



g = E c /A [Uj where, E c , the Thouless energy, is an 
energy scale related to the diffusion time to cross the 
sample and A is the mean level spacing. In a metal, 
E c — hD c i as /L 2 (Ddas — Vpl/d is the classical diffusion 
constant with I the mean free path and vf the velocity of 
the particle) and therefore g cx L d ~ 2 . In an insulator the 
particle is exponentially localized and g oc e~ L ^ where 
£ is the localization length. The OPT is based on the 
following two simple assumptions: a) /3(g) = dl gf 9 ^^ 
is continuous and monotonous, b) The change in the 
conductance with the system size only depends on the 
conductance itself. With this input the OPT predicts 
correctly a MIT for d > 2 characterized by a size in- 
dependent dimensionless conductance g — g c such that 
P(g c ) = and (3(g) > (<) for g > (<) g c . In order 
to make more quantitative statements about the MIT 
it is necessary to understand in detail how the system 
approaches the transition. 

In the case of d = 2 + e (e <C 1) the MIT occurs at 
weak disorder. A rigorous analytical treatment is possi- 
ble by combining diagrammatic perturbation theory and 
field theory techniques [H, 0, [2(|. The results thus 
obtained for critical exponents and critical disorder are 
in agreement with numerical calculations 10]. Later on 
Vollhardt and Wolfle 2l| proposed an extension of this 
theory valid for any d. The idea was to go beyond per- 
turbation theory around the metallic limit by solving a 
self-consistent equation for a renormalized diffusion co- 
efficient. The self-consistent condition was obtained by 
relating ladder diagrams associated with difussons with 
crossed diagrams associated with cooperons. Time rever- 
sal invariance is a required condition for this relation to 
hold. As in the d = 2 + e case, the MIT is induced by the 
growing effect of destructive interference on an otherwise 
metallic state. However a similar self consistent condition 



was also obtained starting from a locator expansion 24 1 . 
Unfortunately some of the results of the Vollhardt- Wolfic 
theory do not agree with recent numerical results: for in- 
stance it is predicted that the upper critical dimension 
for localization is d = 4 and that the critical exponent, 
v, controlling the divergence of the localization length is 
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v = l/(d—2). However numerical results indicate the 
upper critical dimension for localization is d > 6 and 
v « 1.5, 1 in 3d, Ad respectively @, 0- 

From the above discussion it seems clear that except in 
the d = 2 + e limit the MIT cannot be described by any 
perturbation theory around the metallic or the insulator 
side. The ultimate reason for that it is that, according to 
the OPT, the MIT is a fixed point, (3(g c ) = 0, of the mo- 
tion. It is a well known fact that in these situations the 
Hamiltonian may have universal properties completely 
different from the ones observed in the proximity of the 
transition. The above theoretical approaches neither can 
be extrapolated to the physically relevant case of d = 3 
nor predict the dependence with the spatial dimension- 
ality of different parameters describing the transition. 

In order to make progress in this problem a deeper 
knowledge of the dynamics precisely at the transition 
is needed. This paper is a first step in that direction. 
We solve the selfconsistent condition of Vollhardt-Wolfle 
2oL 21 1 but including the spatial dependence of the diffu- 
sion constant predicted by the OPT. With this simple in- 
put we obtain explicit expressions of different parameters 
characterizing the MIT such as critical exponents and the 
slope of the number variance as a function of the spatial 
dimensionality d. The analytical expressions for critical 
exponents and level statistics are in agreement with pre- 
vious numerical calculations 

H H, 0- Throughout the 
paper we assume no interaction, time reversal invariance 
and periodic boundary conditions. The results of this 
letter are thus not applicable to the MIT of the integer 
quantum hall effect [221 ] or in the one occurringin a 2d 
disordered systems with spin-orbit interactions [231 ] . The 
extension of this semiclassical formalism to these cases 
will be published elsewhere [27| . 

One parameter scaling theory and anomalous diffusion 
at the MIT.- In the metallic limit, g — > oo, the dynamics 
of a single particle in a random potential is well described 
by a normal diffusion process. The density of probabil- 
ity P(f, t) of finding a particle, initially at the origin, 
around the position r at time t is described by the so- 
lution of the diffusion equation, P(r,t) 
real space and P(w,g) = 1 



- (2D c!as t) d/2 

iu,+D cla . q * in Fo urier space 
DrinJ,, g oc L d ~ 2 > 1 for 



where \q\ 2 = q 2 . Since (r 2 ) = u clas 
d > 2. These expressions are the starting point to study 
transport properties and level statistics in the metal- 
lic limit or when the localization corrections are small 
[icj |. The MIT for d > 3 occurs in the strong disor- 
der region which is beyond the range of applicability of 
perturbation theory. However the OPT can predict the 
type of motion at the critical point for any dimension: 
at the MIT g = E c /A does not depend on the system 
size. Since A oc 1/L d , the Thoulcss energy must scale as 
E c oc 1/L d . This only can happen if the diffusion at the 
MIT is anomalous with, (r 2m ) ~ t 2m / d j2|| where m is a 
positive integer. This result can likewise be interpreted 



as that the diffusion constant become scale dependent 
D(L) oc 1/L d ~ 2 or, in momentum space, D(q) oc q d ~ 2 . 
The OPT is only capable to predict (r 2 ) but not the dis- 
tribution function P(u),q). However it is precisely this 
function the one needed for the evaluation of critical ex- 
ponents or level statistics. On the other hand it is ev- 
ident that any perturbation theory around the metallic 
or insulator limit will fail if it cannot take into account 
the anomalous diffusion predicted by the OPT. In fact 
anomalous diffusion in low dimensional systems has been 
related to a power-law decay of the eigenstates 0, |26| . 
This strongly suggests [13] that a new basis for the local- 
ization problem given by eigenstates with a power-law 
decay is the starting point for a meaningful perturbation 
theory at the MIT. 

In order to proceed we have to come up with an ex- 
pression for P(oj, q) consistent with the OPT prediction, 
(r 2 ) ~ t 2 l d and that at the same time it can describe 
the dynamics in the proximity of the MIT. The simplest 
alternative is to assume that the classical diffusion coef- 
ficient D c i as gets renormalized to D = D(uj)D(q) with 
D(q) = Doq d ~ 2 . The function D(oj) is given by the solu- 
tion of the self consistent condition of Vollhardt-Wolfle, 
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As a first step to solve this self-consistent equation 
we replace the sum by an integral and write down 
A and D c i as as a function of hp = mvp/% and I, 



D( 
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f/A-: 
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using 7FT ^ 



a d x d-3 
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and noting that the effective dis- 
order strength is controlled by the parameter A = l/kpl 
(the metallic limit corresponds thus with A — > 0), 
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where we have used that the localization length £ 



lim^^o y — D{uo)/iu>. This expression is obtained by 
matching the predictions for the conductivity and den- 
sity response function on the metallic and insulating side 
of the transition 2(J 2l| . 

The third term in Eq.|T]) vanishes (£ — > oo) as we ap- 
proach the MIT. We use this fact to compute the critical 
disorder, A = A c and the critical exponent, s, related to 

the vanishing of the conductivity. The critical disorder 

i 

= { Ji-2d, ) d ~ 1 is obtained by solving A in Eq|T]) with 
lim^^o D{u>) = 0. On the metallic side of the transition 
lim^^o D{uS) = since D(uj) oc o{uj) and the conductiv- 
ity cr(0) oc g c /L = for L — > oo. This result agrees with 



21|. Therefore 



the prediction of Vollhardt and Wolflc 
anomalous diffusion does not affect the behavior of A c 
close to the transition. We note that A c = l/ir > for 
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d — > oo limit. This is consistent with the fact that the 
MIT in the Cayley tree occurs at finite disorder 

The situation is different in the case of the critical ex- 
ponent v related to the divergence of the localization 
length. As we approach the transition from the insu- 
lator side (A > A c ) the localization length diverges as 
£ oc | A — A c |~' y . The conductivity also vanishes on the in- 
sulator side of the transition lim^^o oc iuj 2(| 21|. 
Combining the results for the the metallic and insulator 
sides of the transition it turns out lim^^o D(uf) oc iui = 0. 
Using this fact the critical exponent is, is obtained by 
simply solving Eq.JT]) for £ with w — > 0, 



1 



d-2 



(2) 



A few comments are in order, this expression: a) only 
agrees with the Vollhardt-Wolfle prediction for d ~ 2, b) 
agrees with numerical calculations for any d > 2 dimen- 
sions [7(, c) shows that the upper critical dimension for 
localization is infinity since v > 1/2 for any finite d. This 
is the most important result of the paper. 

Finally we study the critical dimensional conductance 
9c = h L*A • ^ n or der to proceed we have to compute 
D(L). In practical terms this can be carried out by in- 
cluding a lower-cutoff ~ 1/L in the integral over momen- 
tum Eq.JTJ). For A > A c , and uj — > 0, D(uj) — ► and, 

A _Sd_ fl/LJA 2 *- 3 

TT (27T) d JO 



D(L) 



+ q d 



After performing the inte- 



gral, and taking the limit £ — > oo for a fixed L we obtain, 



9c 



n(d- 2)(27r) d 



(3) 



where is the surface of the d sphere. We note that 
g c <C 1 for d ^> 1. This result agrees with pre vious pre- 
dictions based on a self-consistent condition [2l| or simple 
one loop perturbation theory 28] . Thus anomalous diffu- 
sion does not affect the value of g c . However corrections 
to the above g c due to a finite L or £ will be in general 
different from the predictions of Ref. [2l|, Q . 

Level Statistics at the MIT.- We study analytically the 
number variance and the two level correlation function 
(TLCF) at the MIT. Our main result is that the number 
variance is linear S 2 (£) ~ \i with \ w 1 — 2/d. 

We are now interested in the properties of the sys- 
tem precisely at the MIT. Key in our argument is again 
the fact that diffusion is anomalous at the MIT. Our 
starting point is the connected TLCF, iZa(ei,e2) = 
A(p(ei)p(e2)) , (( ) denotes averaging over disorder re- 
alizations and p stands for the spectral density). In 
the metallic limit, g ^> 1 and for s = ei 7 e2 s> g , 
the TLCF is related to P(ui,q ni ) = _^ + j cias ^ by, 

R.2(s) — — ^2$tJ2n - P 2 ( s A,q ni ), where the sum runs 
over all momentum eigenstates q ni . This result is semi- 
classical in the sense that interference corrections rep- 
resented by maximally crossed diagrams are not taken 



into account. As disorder increases diffusion becomes 
slower as a consequence of the growing interferences ef- 
fects. Corrections to the metallic results above are thus 
expected. 

The anomalous diffusion (r 2 ) ~ t 2 / d at the MIT is 
reproduced by simply replacing the standard diffusion 
pole ~ q 2 in P(co,q) by ~ q d . We thus propose that for 
s ^> g c the TLCF at the transition is given by, 



i 



(is + gc\q ni \ d )' 2 



(4) 



where \q ni \ — y X)i=i n f- I n this approach we assume all 
interference effects are included in the renormalization of 
the diffusion coefficient D c [ as — + D q d ~ 2 . Corrections to 
this result are expected clue to to the multifractality of 
the eigenstates. However such corrections cannot modify 
the scale invariance of g c at the MIT. 

We are now ready to compute the number variance, 
= {(N e - (N e )) 2 ) = 2/*(s - £)R 2 (s) with N e the 
number of eigenvalues in an interval of length £ in units 
of the mean level spacing. Carrying out this integral and 
replacing the sum over momenta by an integral, 



E 2 (£) 



1 S d 



dt\t 



id— 1 
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g 2 t 2d 



TT 2 (27T) d J 

Performing the integral and using Eq. ([3]) , 
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A linear number variance with % < 1 is considered a 
signature of a MIT. The origin of this linear behavior was 
predicted heuristically 12| by using OPT and making the 
plausible approximation that eigenvalues interact only if 
their separation is smaller than g c . The value of the 
slope was later estimated to be x — ^J 2 ' EH- We do 
not fully understand the relation between this result and 
ours x = 1 — 2/d. However we note our expression for x 
reproduces correctly the limits d ~ 2 and d — > oo. For 
d > 3 our prediction is around 10 — 15% off the numerical 
value [3], the predictions of Ref. [29| fails for d > 3. 

We now turn to the discussion of the TLCF in the 
region s ~ g c . We aim to examine the heuristic argu- 
ments of Ref. 12| where it was suggested that for s > g c 
there must be a sharp suppression of the spectral cor- 
relations. If localization corrections are negligible the 
use of the supersymmetry method permits an explicit 
evaluation of the nonperturbative part of the TLCF [3(| • 
(a) oc D 2 (s,g) where D(s,g) = EL^oC 1 + ot)" 1 



tdNP 
K 2 



is the spectral determinant associated to the classical dif- 
fusion operator P(uj,q). It is plausible to expect that 
such expression for -R^ p can still be used at the MIT pro- 
vided that the spectral determinant is modified to take 
into account the anomalous diffusion predicted by the 

OPT, namely, \q ni \ = \/£iLi n ? and 9 -> 5c- 



This is 



4 



again a semiclassical approximation, we suggest that all 
the quantum interference effects at the MIT are included 
by an appropriate redefinition of the classical spectral de- 
terminant. The spectral determinant D(s,g c ) can then 
be estimated analytically by exponentiation of the prod- 
uct, replacing sums by integrals and using Eq.©. The 
final result is simply, 

i^ p ( S )oce- 2 ^ d - 2 >/ d (7) 

We note, a) an exponential decay with a similar prefac- 
tor has been observed in numerical calculation 7] , b) the 
sharp decay of the TLCF does not occur for s ~ g c <C 1 
but rather for s ~ 2d/(d — 2)ir 3> g c , c) the conformal 
symmetry predicted by the OPT at the MIT is only con- 
sistent with an exponential or a power-law decay of the 
TLCF, b) the power-law decay of i?2(s) observed in Ref. 
[3l| is related to how the system approach the transition 
rather than to the transition itself. 

Limits of applicab ility .- According to numerical and 
heuristic arguments [32| it is expected that P(q,uj) will 
depend on the multifractal dimension Di rather than d 
for times and distances much shorter than the Heisen- 
berg time and the system size respectively. The proposed 
renormalization of the diffusion coefficient (D(q) ex q d ~ 2 ) 
is in principle restricted to the region q,oj — > 0. However 
we note this region is the only relevant in the calculation 
of the critical exponents. In the evaluation of the number 
variance SQ) other momentum regions may also be rele- 
vant. Therefore our expression of the slope as a function 
of d may get corrections depending on the multifractal 
dimension D2 [29]. These corrections point to the limit 
of applicability of our approach. The OPT is based on 
the scaling of the moments not on the distribution func- 
tion itself. However multifractality is a property of the 
distribution function and consequently beyond the reach 
of the OPT formalism. 

In conclusion, we have put forward an analytical 
approach to the metal-insulator transition based on the 
anomalous diffusion predicted by the OPT at the critical 
point. With this simple input we have shown that the 
upper critical dimension is infinity and found explicit ex- 
pressions for the critical exponents and critical disorder 
as a function of the spatial dimensionality. Moreover we 
have shown that the number variance is asymptotically 
linear with a slope which is a simple increasing function 
of the spatial dimensionality. All our analytical pre- 
dictions are in fair agreement with numerical simulations. 

I thank Denis Basko and Emilio Cuevas for illuminat- 
ing conversations. 
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